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Abstract 

We give a hybrid two stage design which can be useful to estimate the reha- 
bihty of a parallel-series and/or by duality a series-parallel system, when the 
component reliabilities are unknown as well as the total numbers of units al- 
lowed to be tested in each subsystem. When a total sample size is fixed large, 
asymptotic optimality is proved systematically and validated via Monte Carlo 
simulation. 
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1. Introduction 

In reliability engineering two crucial objectives are considered: (1) to 
maximize an estimate of system reliability and (2) to minimize the variance 
of the reliability estimate. Because system designers and users are risk-averse, 
they generally prefer the second objective which leads to a system design with 
a slightly lower reliability estimate but a lower variance of that estimate. In 
the case of parallel-series and/or by duality series-parallel systems just as 
was defined in the variance of the reliability estimate can be lowered 
by allocation of a fixed sample size, while reliability estimate is obtained by 
testing components |2|]. Allocation schemes for estimation with cost jsj lead 
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generally to a discrete optimization problem which can be solved sequentially 
using adaptive designs in a fixed or a Bayesian framework 0, [HI, 0, 0, 
In [8] a reliability sequential schemes (R-SS) was applied successfully to a 
parallel-series systems , when the total number of units to be tested in each 
subsystem was fixed. Recently in joj a two stage design for the same purpose 
was presented and shown to be asymptotically optimal when the subsystems 
sample sizes are fixed and large at the same order of the total sample size of 
the full system. The problem considered in this paper is useful to estimate 
the reliability of a parallel-series and/or by duahty a series-parallel system, 
when the component reliabilities are unknown as well as the total numbers of 
units allowed to be tested in each subsystem. This work improves the results 
in jof by developing a hybrid two stage design to get a dynamic allocation 
between the sample sizes allowed for subsystems and those allowed for their 
components. 

In section |2l definitions and preliminary results are presented accompa- 
nied by the proper two stage design for a parallel subsystem just as was 
defined in and its asymptotic optimality is proved for a fixed and large 
sample size. In section [31 the full parallel-series system is considered and it is 
shown that the variance of its reliability estimate has a lower bound indepen- 
dent of allocation. This leads, in section H] to the main result of this paper 
which lies in the hybrid two stage algorithm and its asymptotic optimality 
for a fixed and large sample size allowed for the full system. In section \5\ 
the results are validated via Monte Carlo simulation and it is shown that 
our algorithm leads asymptotically to the best allocation scheme to reach 
the lower bound of the variance of the reliability estimate. The last section 
is reserved for conclusion and remarks. 



2. Preliminary results 

Consider a system 5* of n subsystems Si, S2, ■ ■ ■ , Sn connected in series, 
each subsystem Sj contains rij components Sij, S2j, ■ ■ ■ , Snjj connected in 
parallel. The system should be referred as parallel-series system (cf. [l|). 
Assume s-independence within and across populations, then the system reli- 
ability is 

n 

R = l[R, (1) 



2 



where 

i?,- = 1 - n (1 - Rij) 

i=l 

is the rehabihty of the parallel subsystem Sj and Rij the reliability of com- 
ponent Sij. An estimator of R is assumed to be the product of sample 
reliabilities 

n 

R = l[R, 

i=i 

where 



Rj^l-\\{1-R,j 



no 

1=1 

and Rij is the sample mean of functioning units in component Sij, 



Ma 

A _ 1=1 
IXij - 



Mij 

Rij is used to estimate Rij where Mij is the sample size and xfj' is the binary 
outcome of the unit / in component S^j. Hence, for each subsystem Sj, one 
must allocate 



i=l 



units such that the estimated reliability of the full system is based on a total 
sample size 

n 

As in the scries case, the variance of the estimated reliability R incurred by 
any allocation scheme can be obtained 

n n 

Var = n (^«^ + - n ^? (2) 



where 



Var 



{a,} 



'l-Ri 



n 



i=l 



(3) 



is given as a function of the allocation numbers Mij and the coefficients of 
variation of Bernoulli populations 



1/R 



ij 



We have found convenient to work with the equivalent expression of 



VarlR 



where 



Rj] 




.-2 



M 



is a sum over all the products of at least two of its arguments. 

The problem is to estimate R when components reliabilities are unknowns 
and a total number of T units must be tested in the full system. The aim 
is to minimize the variance of R. Hence, the problem can be addressed by 
developing allocation schemes to select Mij, the numbers of units to be tested 
in each component i in the subsystem j, under the constraint 



T 



(4) 



=1 1=1 



such that the variance of R is as small as possible. Reliability sequential 
schemes (R-SS) exist for the series, parallel or parallel-series configurations 
when the sample sizes Tj of the subsystems are fixed. Therefore, one can 
fully optimize the variance of R just by applying the (R-SS) to find the best 
partition Ti,T2, ...,T„ of T. Unfortunately, a full sequential design can not 
be used in practice for large systems since the number of operations will 
growth dramatically. For this reason, we reasonably propose a hybrid two 
stage design which is shown to be asymptotically optimal when T is large. 
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2.1. Lower bound for the variance of the estimated reliability of the parallel 
subsystem Sj 

For the asymptotic optimization of the variance of the estimated rehabih- 
ties, we make use of the well-known Lagrange's identity which can be written 
in the form: 

Let Oi > 0, Ni > 0, for i = 1, A; and = A^i H h A^fc, then the 

following identity holds. 



1=1 

Proposition 1. Denote by 
then 

Var \Rj\>Q 



i=l j=i+l 



j=i 



Proof. The proof is a direct consequence of the previous identity 
deed 

2 

/ '"J \ 

Var 



ra,— 1 n, 



i=l k=i+l 
-2 -2 -2 



(5) 



(6) 



In- 



(7) 



2.2. The two stage design for the parallel subsystem Sj 

Following the expansion ([71) and since F contains second order terms (see 
later), one gives interest to the numbers Mij which minimize the expression 



rii — l Ui 



n2 



Tri J2 E ~ 



i=l k=i+l 



M.jMkj 



Thus Mij must verify for i 
which imphes that 



l,...,nj 



T - 



(8) 



k=l 

If one assumes that Tj is fixed then a proper two stage scheme can be used 
to determine Mij, just as was defined in j9[, as follows: 
Choose L function of Tj such that: 

(i) Lj must be large if Tj is large, 



(iii) lim ^ 



0. 



One can take for example Lj = [^/Tjj , where [.] denotes the integer part. 

Stage 1. Sample Lj units from each component i in the subsystem j, esti- 
mate Cij by its maximum likelihood estimator (M.L.E) 



\ 1=1 

and define the predictor, according to ([8]) 



T- 



fc=i 



l,...,nj - 1 



Stage 2. Sample Tj—rijLj units for which Mij—Lj are units from component 
i in the subsystem j where M^j is the corrector of defined by 



M, 



■^i3 



max i Lj , Mij > , z 



fe=i 



l,...,nj - 1, 
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Theorem 1. Choosing the Mij according to the previous two stage sampling 
scheme, one obtains 



lim Tj Var \Rj \ -Q 







Proof. From relation ([7j), one can write 



i=l k=i+l 



T, (Var - Q,) = (1 - Rjf E 



(9) 



When Tj is large enough, condition (iii) gives Mij = Mij for i = 1, rij — 1. 
So, the strong law of large numbers with the integer part properties give, 
when Tj — )■ oo, 

Mkj Cij ' 
for i = 1, ...,nj. Hence, 

-^2 



MijC,} - MkjC,/) _ M,, / _i Mkj _i 



and on the other hand 



c"> ^0, as Tj oo, (10) 



(11) 



which achieves the proof. 



3. Lower bound for the variance of the estimated reliability of S 

We consider now the full system S. From expression one can write 



Var\R[ =R 



" / Var ( R 

n 



R] 



1-1 



The following theorem gives a lower bound for the variance of R. 
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Theorem 2. Denote by 
Q = T-^R^ 

then 



El — Rj 



Var > Q 

Proof. Expanding the right hand side of and using ([T]), one obtain 



VarlR \ =R 



E 



Var R 



+ F 



Var Ri Var R 



Rl 



i?2 



which gives with the help of theorem ([T]) 



j = l 3 

This last expression has the form 

n 



which can be expanded, thanks to identity ([5]), as follows 



2rr-l 



R^T 



El — Rj 



R 



;=1 ^ \k=l 



i=l j=i+l 

and as a consequence 



3 

n-l n \J^i Ri Z-^ Hi 



k=l 



k=l 



TT- 



Var > T-^R^ 
which achieves the proof. 



.j=l ^ \k=l 



Q, 
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4. The hybrid two stage design for the fuU system S 

Similarly to the case of a subsystem an from expressions (IT^ and (fT5|l . 
one gives interest to the numbers Tj which minimize the quantity 

2 



n— 1 n 

EE 

i=l j=i+l 



k=l 



k=l 



and obtains the asymptotic optimality criteria 



1^ V 



k=l 



rp Ui 



k=l 



for all i, j G {1, 2, n}, which gives the rule 



1-_R. 



T. = T 



fc=l 



Z^ Wfc 



(14) 



We can now implement a hybrid two stage design for the determination of 



the numbers Tj as well as Mij as follows: 



Stage 1 choose L 



Vt 



: one applies the two stage scheme given in sub- 
section (12.11) for each subsystem Sj with Tj = L and Lj = [\/Tjj . 
Next, obtain the predictor, according to the rule flT4|) . 



T- 



l-Ri 



fc=l i=l 



n — 1. 



Stage 2 define the corrector 

Tj = max|L,T;- [> , j = 1, . . . ,n - 1, 

T 



n-l 



and take back the two stage scheme for each subsystem Sj to calculate 
Mij with the sample size equals Tj. 

Now, the main result of this paper is given by the following theorem. 

Theorem 3. Choosing the Tj and Mij according to the hybrid two stage 
design, one obtains 

\mv T {Var^R^-Q^ = 

where Q is defined in theorem 
Proof. The relation (jH]) implies that 



Var{RA=Q,+Tr^{l-R,f t 

I — 1 L—'i-l 



MijMkj 



+ {l-RjfF 



i=l k=i+l 



As a consequence of the hybrid two stage design and the strong law of large 
numbers, T /Tj and Tj/Mij remain bounded for all z, j as T — t- oo. It follows 
that, as T — i- oo, 

and 

i=i k=i+i '^^ 
thanks to ([lOD and (HH). Thus, 

Var = Qj + o (T^^) , as T -> oo 

which implies that 




Var-!^Rj 



+ 1 



n(|.i.or-)) 



i=i ^ J 
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As a consequence 



Urn T {VariR} -Q] = R-" lim T. 

r-i>oo V I J / T-s>oo 



.7 = 1 ^ 



n 



+ 1 -1-g 



Now, expanding the product within the hmit and applying identity dSJ, after 
having replaced Qj by its expression ([H]), one obtains 



• 1 v^' 



where 



^2 + ^ I i?2'-'i?2 



i?2 



,2 



<5l Qn 



n— 1 n 



A 



B 



7^. 



Rk 



1 = 1 



■n-k \ 



1=1 



1=1 k=i+l 



Once more, the hybrid two stage allocation scheme and the strong law of 
large numbers provide 



and 



lim T.A = 



lim T.B = 



which achieves the proof. 



5. Monte Carlo simulation 

Let us remark first that the lower bound Q is a first order approximation 
of the optimal variance of the reliability estimate under the constraint (jl]) 
when T is large. 

In the first experiment, as in figure ([1]), we consider a simple parallel- 
series system of two subsystems each one, with varying reliabilities and a 



11 



fixed sample size T = 20. For each situation A,B,C and D and for eacfi 



Vt 



to T 



Vt 



partition sample size {Ti, T — Ti} wliere Ti varies from 

we have apphed the proper two stage design for each parallel subsystem 

and reported in a bar diagram Var (jij as a function of Ti, see figure ([3]). 

On the other hand, in table ([T]), we have reported the expected value of 
Ti = Mil + M21 given by the hybrid two stage design. As expected, our 
scheme gives the best allocation for each situation. 

The second experiment deals with a non trivial parallel-series system 
just as in jof, where subsystems are composed, respectively, of 2,3,4 and 5 
components, see figure ([2]). The partition total numbers Tj to test in each 
subsystem are evaluated systematically by the hybrid two stage design while 
their sum T is incremented from 100 to 10000 by step of 100. The figure 

@ shows the rate of the excess of variance T (^ar (jij — at logarithmic 

scale as a function of the sample size T. The asymptotic optimality of the 
hybrid scheme is validated. 



Figure 1: A simple parallel-series system of two subsystems with two components each 
one 



6. Conclusion 



The proof of the first order asymptotic optimality for the proper two stage 
design for a parallel subsystem as well as for the hybrid two stage design for 
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R,,=0.9 



R,,=0.1 



R,.=o.i 



"13=0-3 



"23=0-1 



"33=0-8 



"34=0-1 



44 



"43=0-8 



54 



Figure 2: A non trivial parallel-series system. 



System 


— R21 — R12 — R22 




A 


0.1-0.11-0.9-0.99 


16 


B 


0.5-0.55-0.51-0.6 


11 


C 


0.9-0.99-0.1-0.11 


4 


D 


0.2 - 0.4 - 0.6 - 0.3 


12 



Table 1: Expected value of Ti = Mn + M21 given by the hybrid t-wo stage design 
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Figure 3: Bar diagram Var (^R^ as a function of Ti for each case A,B,C and D : "shows 
the minimum of Var (^Rj 
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Figure 4: Asymptotic optimality of the hybrid two stage design : the speed of the excess 
of variance T ( Var (r) — Q] at logarithmic scale as a function of the sample size T 



the full system has been obtained mainly through the foUo'wing steps 

• an adequate writing of the variance of the reliability estimate, 

• a lo'wer bound for this variance, independent of allocation, 

• the allocation defined by the hybrid sampling scheme and the strong 
law of large numbers. 

With a straightforward but tedious adaptation, the above study can be 
namely extended to deal with complex systems involving a multi-criteria 
optimization problem under a set of constraints such as risk, system weight, 
cost, performance and others, in a fixed or in a Bayesian framework. 
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